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Abstract. 

This article describes a simple heat model of the vacuum arc remelting (VAR) 
process that includes solution of the nonlinear heat conductivity equation with the 
nonlinear boundary conditions which are typical for VAR process. The finite-difference 
analogue of the model equations was obtained through the finite volume method. To 
check the efficiency of the simplified model that does not include magnetohydrodynamic 
phenomena in the liquid metal pool, the comparison has been made of the numerical 
calculation of the metal pool depth when melting the Russian titanium alloy VT3-1 
with the results of radiographical tests. It was established that the model adequately 
describes the test data for various melting modes (ingot diameter and current strength). 



1. Introduction. 

Today a lot of attention is given to automation of VAR furnaces and the focus is 
on development of the arc clearance control flowcharts, feedback controllers, program 
management systems for furnace electric modes pQ. Together with the aforementioned 
issues, the development of the VAR automatic control systems requires a reliable 
theoretical description of the ingot formation process that underlies the automatic 
control system algorithms. The significance of through and reliable theoretical models 
of the ingot process formation is determined by the fact that the immediate control 
of the ingot parameters during solidification of titanium alloy (in particular, the depth 
of the metal pool, width of the mushy zone, etc.) using the measuring instruments, 
appears to be a problem so far. Development of the theoretical ideas about formation 
of titanium alloy ingots is rather important today since there are cases when different 
ingot-related defects are observed on billets, bars and eventually on critical parts. Such 
defects include /3-flecks [2], zonal segregation [3], tree-like structure [I] and dark/light 
spots [3] that may be observed in the transverse section under the macro examination 
of mill products in different titanium alloys. Very often such defects are related to 
the final remelting stage therefore the theoretical description of the processes that take 
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place during solidification of the VAR melted ingots, as well as understanding of the 
impact of various melting parameters on ingot metallurgical quality, is undoubtedly the 
paramount objective for creation of the fully automated VAR control system. 

Many researchers studied the heat processes in ingots during VAR [HI IB], 
electroslag remelt |Hl E] and continuous casting [TU] . Today, in view of development 
of advanced computers it is getting important to write science-based VAR simulation 
programs for their use in industrial environment in order to develop new remelting 
modes for complex titanium alloys. 

Taking into account aforementioned facts we reviewed in this article the 
mathematical model of the vacuum arc remelting and compared the theoretical 
calculation of the metal pool depth with the radiographical test results obtained for 
Russian titanium alloy VT3-1 (Ti-6.5Al-2.5Mo-l.5Cr-0.5Fe-0.3Si) p]. Taking into 
consideration the non-trivial complexity and essential non-linearity of the task, it should 
be noted that firstly there is no exact analytical solution of the task and secondly today 
there are a lot of different methods of approximation of the boundary conditions and 
building of the difference schemes for the heat conductivity equation. It is obvious 
that in such situation the experimental check of the computational solutions obtained 
through different methods of approximation is getting important. Therefore, in this 
article we reviewed both aspects of the problem: making the mathematical model and 
comparison of this model with the experimental results. 

2. Vacuum arc remelting process. 

We will describe the heat flows during the vacuum arc remelting through the heat 
conductivity equation that includes release of latent heat within the interval of the alloy 
solidification. In order not to involve the complex interface conditions on the phase 
boundary, we reviewed the so-called method of apparent heat capacity [TH] that allows 
to include crystallization rate (i.e. heat source) into the system heat capacity as a 
complementary additive term [5]. Such problem statement allows to make calculation 
in one domain without "distinguishing" the liquid, solid and mushy zone. 

Also, from this point on we considered that density, heat capacity, heat conductivity 
and heat-transfer coefficient depend on the temperature in an arbitrary way [T2"j . 

Heat transfer equation for VAR ingot can be written in a form 



where T is temperature, r is time, C{T) is heat capacity, A(T) is heat conductivity, 
p(T) is density, r and z are radial and axial coordinates. 

The boundary may be divided into several areas shown in Figure H Such division 
is based on physical processes that take place on the ingot surface under VAR. We will 
consider those processes and the relevant boundary conditions. 




(1) 
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Boundary AB. On the ingot axis we have the simple symmetry condition: 
dT 

*=°- (2) 

Boundary BC . The border BC corresponds to a bath mirror surface which are being 
under an end face of the consumable electrode. On this border we shall set temperature. 
In the elementary kind the temperature of a surface of a bath is defined by an overheat 
above alloy liquidus temperature 

T = T L + AT, (3) 

where Tl is liquidus temperature, AT is overheat. 



Boundary CD. The boundary CD corresponds to a ring gap. Here we also set a 
boundary condition of the I-st kind 

T = T L + ° cr ~ 2 n r -AT, (4) 

where D cr is crucible diameter, D e i is electrod diameter. 



Boundary DE. The given boundary corresponds to a zone of a contact belt for which 
experimental values of a specific thermal flux q ^3] are known, therefore we write down 
a boundary condition of Il-nd kind 

dT 

-A(2% = g. (5) 



Boundary EF. The site EF corresponds to a zone of forming ingot on which has 
already partially or completely occured separation of solid ingot from a crucible surface 
due to passage of process volumetric shrinkage. 

Thus three possible mechanisms of heat removal can be observed: 

• Heat removal by radiation through the formed gap. 

• Heat removal by contact way in places of contact of a surface of an ingot with an 
internal wall of copper crucible. 

• The heat transfer by convection through a gas phase in a gap. 

Three specified mechanisms, it is possible to write down the boundary condition 
considering all in the form of 



- ACT)— = rj (e<r (T 4 - T 4 ) + - (T - T cr )j + (1 - n)a (T - T cr ) , (6) 

where n - factor describing a share of contact of a surface of an ingot with a crucible, 
e - emissivity, T cr is temperature of inner surface of the crucible, a is heat transfer 
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Figure 1. Vacuum arc remelting process. 



coefficient in the contact zone, k is heat conductivity of gas in gap, d is gap. The 
parameters r], e, T cr , a, k d are depending on temperature. For simplification of a 
record we shall indicate a following designation 

At CO = ( ea - _ ^ + - j + (1 - T]cr)0Ccr, (7) 

and © is rewritted in a more compact form 

-A(T)^=/3 cr (T)(T-T cr ). (8) 

Thus, complex process of the heat transfer on border EF is described with the nonlinear 
law JHJ. 

Boundary FA. On a site FA boundary condition is similar to (JHJ), therefore we shall 
write down: 

dT 

-X(T)— = (3 bot (T)(T-T bot ), (9) 
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where 

MT) = Vbot (ea T ~^ bot r, bot )a bot . (10) 

\ J- J- bot a J 

Besides it is necessary for us to set the initial condition for the equation (JTJ. During 
vacuum arc remelting melt an ingot occurs, i.e. the border BCD moves upwards 
with the speed defined by remelting conditions. Melting rate of an ingot is defined 
by the mass speed of fusion depending on force of a current, an electrode and a filled 
crystallizer diameters, and also on some other parameters. We shall accept in the 
beginning of remelting is already available ingot of some small height with homogeneous 
distribution of temperature on all volume is considered to be as the initial condition. 
Initial temperature it is accepted equal Tl + AT. For the subsequent moments of 
time the initial condition is defined by distribution of temperature during the previous 
moment of time with additional small layer of an ingot which temperature is Tl + AT. 

In the following sections we shall consider discretization of the given problem on a 
rectangular grid, including linearization of a system, and decision method of the received 
the linear equations system. 



3. Discretization of heat conduction equation. 

We will consider the most common case when A(T) and p(T)C(T) are piecewise 
continuous functions. In this case it is convenient to use the integral identity method 
to the finite-difference scheme construct ^2] . We introduce the heat fluxes along r and 
z directions: 



dT 

Q(r,z) = X(r,z) — , 



P(r,z) 



A(r, z 



dT 

dz 



(12) 



Taking into account these notations, Eq. (JTJ) is rewrited to 
p(T)C(T)^ = ll(rQ(r, Z )) + j-(P(r,z)). 

Integrating the Eq. (|T3~|) over the cylindric volume for r G [^-1/2,^1+1/2] and z G 
[ z j-i/2, z j+i/2\ ( see Figure El) we can write heat balance equation can be put as follows 



(13) 



r i-l/2 z j-l/2 



rdr I dzp(r, z)C(r, z) 



dT(r, z) 
dr 



r i+l/2 z j+l/2 



2=3 + 1/2 



^+1/2 



+ 



n+i/2Q{n+i/2, z) - r i -. 1 /2Q(r i - 1 /2, z) 
P(r, z j+1/2 ) - P{r, Zj-1/2) rdr. 



dz 



(14) 



-1/2 
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This equation is exact. Now we need to evaluate the integrals, which contains Eq. (|14|) . 
Because we have (JTTJ) than 

Q[r ' Z) dr = T(r l+1 )-T(r i ). (15) 



A(r, z 
On the other hand 



rM \ ri+1 j 

Q(r, z) , , f dr . . 

Y^-rdr « Q(r i+1/2 , z) / -— — . 16 

Air, 2) J Air, z) 

Thus we can write following expressions within the fluxes for internal region of ingot 
Q[r i+ x/ 2 ,z)tt , (17) 



/ 


dr 
\(r,z) 


r(r f , ^) - 






dr 


/ 

n-i 


\(r,z) 


T ( r ; Zj+l) 




Zj + l 

I 

Z 3 


dr 
A(ry») 


T(r, zj) - 


- T(r , Zj^x) 



p(r, z j+1/2 ) « v ' c;; — (19) 



f(r, ^i-1/2) - V ' % lAil^Z . (20 ) 

r cir 



J A(r, 2 ) 
z,_i 



Let these expressions be substituted in Eq. (JT3j) then 



' ' 1 Zj ~ 1/2 dT(r z) 

rdr J dzp(r, z)C(r, z) — 7^ — 

i+l/2 z j + l/2 



Z 3+l/2 , 2 J + l/2 



r i+l/2 - [ — x r i-l/2 [Ti,j - Ti_ij] I 

J c rlr J 



r dr J f dr 

Z i-V 2 J \fr~z) Z i-V* J \(r,z) 

ri ri-! 

r «+l/2 r i+l/2 

+ Pm' + i-T m ] / 1 JLJ— + [T U -T id - 1 ] j (21) 

J r dz J r dz 

ri -v 2 J A(7^y r *-va J A(7^y 

Zj z j — l 

Making use the simplest approximation for the integral in the left side of equation we 
obtain 
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\ 'i+1/2 ' i-l/2 n 
Zj-l/2 - Zj+l/2) 7, Pi,j U i,j 



T 

z j + l/2 z j + l/2 

dz r „ „ , r dz 



n+1/2 [Ti+i,j - Tij] I — rj_i/ 2 [Tij - Ti-ij] I 

J r dr J f dr 



r dr 

n y > ; ri-i v ' 

n+i/2 n+1/2 
rm m i f fdr rm m , f rdr , . 

[Tij+i-Tij] / — + [Ty-r w _ a ] / - . (22) 



r dz J r dz 



Z j Zj _ 1 



Moreover the integrals containing heat conductivity is simplified up to: 

Zj+l/2 



\+l/2,j, 

r dr i i+1 — ' i 

Zj-l/2 
Zj+l/2 



dz 


_ Zj+l/2 - Zj-1/2 


n+i 

r dr 
I A(r,«) 


n+i - u 




_ Zj+l/2 - Zj-1/2 


r J , A(r,z) 




rdr 


_ 'i+1/2 'i-l/2 


r dz 
J A(r,z) 


2(Zj+l - 


rdr 


/>-• 2 ,yi 2 

^ ' i+1/2 ' i-l/2 


z j 

r dz 
J A(r,*) 


2 (^i - Zj-l) 



i+l/2 

nil- I" , ,, — /T , 

AjJ+l/2) 

j , ,/- _i. , ,_(_ ] •- / / 

n-1/2 

r »+l/2 

;■(•//■ , .-> — l'~ , .-, 

Aij-1/2- 

r dz ~~ *i-l) 

-1/2 



After that the Eq. f!22l) can be rewritten as 
Ta ■ — T- ■ 



r 



r 



i+1/2 



2 2 
''i+1 r « 'i+1/2 'i-l/2 



/t-i T-i \ \ r i-l/2 g 

l-'ij ~~ J *-ljJ ' A i-l/2j • " ~2 ~2 



+ (+i,j+l — ^i,j) ' ^i,j+l/2 



iy-i iy /yi ^ yi — 

'« '*-l 'i+1/2 'i-l/2 
1 

(Z j+ i - Zj) ■ (zj+l/2 - Zj-l/2) 



- (Tij - Tij_i) ■ Xij-1/2 ■ ~ —7 r- (23) 

\ z i ~ z j-i) " V z i+i/2 ~ z j-i/2) 



We define (for the simplest way) 



_ n + r i+1 _ r 4 _i + r< 
"i+1/2 — ^ > r»_i/2 — : 
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Zj + Zj+i z j — l i z j 
Zj+l/2 = ~ , Z j-l/2 = ^ ' 

Thus for the rectangular grid Eq. (|23|) can be written in the common form as 



T- ■ — Ta ■ 

± i,3 A h3 

Pi,3 U i,j 



r 



_ r i+ i + r j 4A w / 2 j , _ , 

Ti+i-ri ( n+1 + nf - (n - ri-if 

n + rj-i 4Ai_i/ 2 j \t -T 1 

n-n^ ( rm + ri ) 2 -( r .- ri _ 1 ) 2 

2 -^,j + l/2 rrp rp 1 

_ z j) ' l^i+i _ 



2Ajj_i/2 



For constant space steps (h r and /i 2 ) we have 



Fu - . (24) 



T- —T 



r 1,3^1,3 



^ • (l + ±) \T M , - T,J - .U-- 



r 



hi V 1 + 2i ) [Tl+1 " ~ Tl >> ] ~ ' I 1 " 2i) [TiJ " ^" 1JJ 



• [T W ' T id ] - ^ • [T M - Tij-i] . (25) 
To evaluate heat conductivity at the semi - integer domian points we use approximation 

2AjjAj + ij 2Aj_ijAjj 



*i+l/2,j - T — : , A i-l/2,j 



_ 2\jj\ iij+ i _ 2X itj _ 1 X iij 

M,j+l/2 - T — T , \j-l/2 - T — T • (/Oj 

Equation ()25j) joins the temperature Tjj with Ti_ij, T i+1 j, and than Eq. 

(125(1 has 5-point space scheme and explicit scheme in time. We will use 2-cycle scheme 
|12j . Let we rewrite Eq. ()25)1 as follows 



T { ■ — y; ■ 

Pi,fiij — '■ — — — — — A r Tjj- + A z Tjj, (27) 

where operators A r 2 are defined by Eq. (|25j). The main idea of the 2-cycle factorization 
scheme is making of the following steps (A; is time step index) ^2] 

rpk-l/2 _ rpk-1 rpk-1/2 _|_ j,fc_l 
Pi,j C i,j— — = K— n ~ ( 28 ) 
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rpk rpk—1/2 rpk _i_ rpk — 1/2 

P,/', " J"' = A, iJ \ iJ (29) 

rpk + 1/2 rpk rpk + 1/2 rpk 

,>,./',, " T " = \ " , " (so) 

rpk + l _ rpk + 1/2 rpk + 1 , rpk + 1/2 

" r j = X " , " (31) 

The cycle of evaluation is namely to solve equations (J2HJ) - d3H) - Let we rewite these 
equations in a more useful form 



[BiA Tl\f - [A itj + By + 2p i:j C,j/r} T^ 1/2 + [Aij] T^j 2 

\ B i,j\ Tttj - [Aj + By - 2/HjdjM + [Ay] T££\ (32) 



[Ay] I* _ x - [Ay + T M + 2 P yCy/r] T* + [Ty] J* +1 

= - ([Ay] T'jll 2 - [Ay- + ry - 2pyCy./r] ^ 1/2 + [Ty] Tt; + f) (33) 



[By] if^f 



[Ay + Bij + 2pijCij/r] Tij ^ + [Ay] Ti+ij 



- ([By] Ity - [Aj + By - 2pyC y /r] + [Aj] 



(34) 



[A,,] rfit 1 ! - [Ay + r w + ■>,>,,(■„: t\ T^ 1 + [Ty] Tfift 

'[Ay] - [Ay- + ry - 2 P yCy/r] I# 1/2 + [Ty] ) (35) 



where we denoted 



N Ay-1/2 r _ Ay+i/g 



The solving of these equations is strightforward, because equation matrixes are three- 
diagonal, then these equations have 3-point space scheme. It is a very important note 
for our method, because we will have possibility to decrease CPU time. 

It might to prove 2-cycle method is absolutly time step stability and it has 
second-order in r 2 like Crank-Nicholson scheme and second-order in space. Moreover 
defferential operators A rz in the 2-cycle method can be noncommutating ones. 
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4. Discretization of boundary conditions. 

To solve the equations (j3*2j) -(j3*H j) it is necessary to define values of temperature on 
border or to set some relationships connecting temperature on border of a body with 
temperature of an environment. Boundary conditions in a continuous limit are certain 
by equations (J2J) - (jUJ). We shall construct their discrete finite difference analogues on the 
uniform grid, using a method of finite volume [Hj. For this purpose we shall integrate 
the equation ()13j) on some volume G, adjoining to border of an ingot and on time within 
the limits of [t, t + r]. Such volumes are shown on Fig. |21 



j+1 — 

I 

j-i — 



i-l i i+1 



N 



Figure 2. The volumes for evaluating of integral balances in the boundary conditions. 

Integral balance of heat can be written as 

t+ f T r dT T f (l d dP\ 

1 1 omcm^™ = / / + ^ dvdt, (38) 

t G t G ' 

where heat fluxes are defined by Eqs. (fTTj) and (JT2J). 

After integration we shall receive the relationship connecting a thermal flux on 
border of a body, set by a corresponding boundary condition, and thermal fluxes in 
nearby points for which some approximations are required. Integration on time is led 
as follows: the left part of the Eq. (J38|) we approximate as follows 

j j p(T)C(T)—dVdt~T-T, (39) 

t G 
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and the right part of Eq. (J38)) we shall calculate in previous the moment of time. Now 
it is easy to calculate finite difference counterparts of boundary conditions. 

4-1. Boundary AB. 

Let area G = [0,h r /2] x [(J — l/2)h z ,(j + l/2)h z ] then the general expression for a 
boundary condition on a site AB after integration can be written down in the form of 

h 2 h z r - 1 rh z 

— — pC ,j T 0>j - Tqj = —\l/ 2 ,j [Tij - Tqj] 



8 



rh 2 

+ g^ 1 (*0j+l/2 Poj+i - Tqj] ~ Aoj-1/2 [T ,j - Tbj-j) . (40) 



4-2. Boundary BC . 

On a site BC there is no necessity to use a method of finite volume as in this case the 
field of temperatures is set directly, therefore 

T m = T L + AT. (41) 
4.3. Boundary CD. 

On a site CD we set the temperature, as well as on BC 

D cr — 2ih 
D„-D, 



,w =T L + n cr 77 AT. (42) 

cr 



4-4- Boundary DE. 

Let G = [(N — i)h r , Nh r ) x [(j — l/2)h z , (j + l/2)h z ] then the general expression for 
a boundary condition on a site DE after integration it is possible to write down in the 
form of 



h;h z 



(4iV - l)pC NJ f NJ - T NJ = -rh z Nh r q - rh z (N - -)X N -i/ 2 j [T N ,j - Tv-i,j] 
rh 2 

+ gfo~(^N — 1) (Ajvj+i/2 [^ivj+i — Tivj] — Ajvj_i/2 [Tjvj — Tjvj-iJJ • (43) 
^.5. Boundary EE. 

Let G = [(AT — l)/i r , A^/i r ] x [(j — l/2)h z , (j + 1/2) /i^] then the general expression for a 
boundary condition on a site -EF after integration it is possible to write down area in 
the form of 
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= - Th z Nh r (3% tj [T NJ - T cr ] + rh z (N - -)Ajv-i/ 2 j [T Nij - T N _ hj ] 
rh 2 

+ ^(4iV - 1) (\ N , j+ i /2 [T Nd+ i - T NJ ] - Xnj-i/2 [T Nd - r^-i]) (44) 
4-6. Boundary AF. 

Let G — [(i — l/2)h r , (i + l/2)h r ] x [0, h z /2] then the general expression for a boundary 
condition on a site FA after integration it is possible to write down area in the form of 



ih 2 r h z 

n P^ifl 



irtit 



(•\,i/2 [TiA — T i)0 ] — h z (3 i [T ifi — T] 



h z 

= rh z ((i + A i+ i/2,o [7i+i,o - TiA ~ ~ g) ^-V^o P^o - ^l-i,o]^) (45) 

Further we use approximations for heat conductivity factors in semi-integer points, 
similar to Eq. 



5. The model parameters for VT3-1 alloy. 

As it is possible to see from expressions for boundary conditions on a mirror of a liquid 
bath (BC and CD) we set distribution of temperature which is defined by alloy liquidus 
temperature and an overheat depending on parameters of VAR. The overheat of the melt 
A can be described by formula [Sj: 

AT(J,A„) = 400e- 12 ^, (46) 

where J is arc current, kA; D in - ingot diameter, m. In a zone DE the thermal flux 
q to which measurements for alloy VT3-1 work [H] is devoted is set. The parameter n 
defines the relative contribution of each of mechanisms of heat removal: at rj — 1 the 
heat-conducting path goes only by radiation, and at rj = only by convection. For EF 
we accepted T]ef = 1.0, and for FA - t/fa = 0.5. 

C(T), including rate of release of latent heat in an interval of temperatures between 
liquidus T L and solidus T s , it is possible to present an effective specific thermal capacity 
in the form of jlUj : 

C L (T) for T>T L 

C(T) = { g{T)C s {T) + [1 - g(T)\C L {T) - L d -^ for T s <T <T L (47) 
C S (T) for T<T S 
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where Cl and Cs are specific heat capacities for liquid and solid phases; L is latent heat 
of fusion, g(T) is solid fraction. In case of binary alloy for g(T) we can obtain simple 
expression (under lever rule suggestion) 

g(T) = Tm ~ Ts . ILZZ. ( 48 ) 
y{ ' T L -T S T m -T y } 

For temperature derivative of g(T) we can obtain 

dg{T) T m -T L T m -T s 
dT ~ Tl — Ts (T m — T) 2 ' 1 > 

where T m is fusion temperature of a solvent, Tl is liquidus temperature, Ts is solidus 
temperature. Parameters of model for alloy VT3-1 have been chosen by the following 
T m = 1668 °C, T s = 1550 °C, T L = 1620 °C, L = 355000 J/kg, T cr = 70 °C, 
Tbot = 70 °C, abot = 300 W/m 2 K. Some data can be found in [T3]. Position of a point 
E was defined as in work [Bj. 

For the account of increase in heat conductivity due to fluid flow we have increased 
heat conductivity of an alloy in a liquid phase 0. Temperature dependence of the 
resulted ingot surface emissivity degree and an internal surface of crucible has been 
chosen in the form of square-law dependence on temperature . 



6. Simulation results and radiographical experiments. 

Using the mathematical model described above, we have led solidification modelling of 
ingots from the alloy VT3-1. Modes of remeltings are specified in the Tabled 

Some works [3 E] describe the experimental research of the metal pool depth and 
profile through fixation with radioactive isotopes. Below we compare the simulation 
results with the experimental data on liquid metal pool depth at different points of time. 
In order to adequately describe the experimental data we have carried out "adjusting" 
of the model. The adjusting provided for rather exact match of the calculated and 
experimental metal pool depth, as well as liquid metal pool profiles at different ingot 
height. Due to the fact that the most extensive experimental data are available for 
750 mm dia ingot melted at 37 kA, we carried out "adjusting" for the above-mentioned 
melting mode. The main variable parameters were - heat conductivity in the liquid 
phase (that was considered not dependent on temperature) and the surface emissivity 
factor depending on temperature. 



Table 1. Remelting regimes. 



Ingot diameter, mm 


Arc current, kA 


750 


37 


570 


25 


435 


15 
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6.1. 0750 mm ingot. Arc current J = 37 kA. 

One of the most important parameters assessed during analysis of one or another VAR 
mode for titanium alloys (as well as for nickel-base, iron-base and zirconium base alloys 
[3 13 El EH]) is the depth of the liquid metal pool. Figure El shows the theoretical 
metal pool depth with experimentally obtained depth (through fixation with tungsten 
radioactive isotopes) depending on the height of the melted ingot. As is obviously, 
the theory describes the metal pool depth behavior satisfactorily. During melting the 
quasi-steady state was achieved at the metal pool depth having been steady. 

3000 
2500 
2000 

E 

E 1500 
x~ 

1000 
500 


500 1000 1500 2000 2500 3000 3500 4000 
z, mm 

Figure 3. Experimental and theoretical pool depth versus the height of the melted 
ingot. 0750 mm ingot, arc current J = 37 kA. Points - experiment, curve - theory. 

The description of the liquid metal pool only may be misleading as the pool depth 
does not actually reflect the volume of the liquid metal at the particular melt point. In 
order to assess feasibility of satisfactory simulation of the liquid metal pool volume using 
the model, we calculated the liquid metal pool profiles at the points of time recorded 
on the experimental radiogram. This shows that the model adequately simulates the 
profiles of the liquid metal pool during the whole melt process. All this together suggests 
that on the basis of the simulation approach in question it is feasible to simulate some 
most important parameters of metal solidification under VAR (width of mushy zone, 
temperature gradient, isotherm travel rate). 

Figure |U shows the calculated profiles of the liquid metal pool during first 140 
minutes of melting (with 20 minute interval) to evaluate the linear solid-melt interface 
travel rate in different ingot zones. The comparison of the simulated profiles with the 
experimentally obtained metal pool profiles shows that the simulated metal pool volumes 
do not deviate from the experimental ones by more than 15%. 
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Figure 4. Calculated profiles of the liquid metal pool during first 140 minutes of 
melting (with 20 minute interval). 

6.2. 0570 mm ingot. Arc current J = 25 kA. 

Figure El shows the theoretical curves and experimental points for 570 mm ingot melted 
at 25 kA. In this case we compared the simulated data with the experimental data 
according to location of the pool open surface (i.e. melt mirror) and pool bottom 
coordinate versus the melting time. It can be seen from Figure El that the quasi-steady 
state was not achieved. 

6.3. 0435 mm ingot. Arc current J = 15 kA. 

Figure El shows the theoretical curves and experimental points for 435 mm ingot melted 
at 15 kA. In this above we compared the simulated data with the experimental 

data according to location of the pool open surface (i.e. melt mirror) and pool bottom 
coordinate versus the melting time. It can be seen from Figure El that the quasi-steady 
state was achieved. 

Therefore, a conclusion can be made that under selected parameters for Vt3-1 alloy, 
the simulation model can describe the main behavior tendencies of such VAR parameters 
as liquid metal pool depth and profile at different melting points. It is very important 
to note that the model that was "adjusted" once, gives satisfactory description of the 
experimental data beyond the "adjusting values". 
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Figure 5. Location of the pool open surface (upper curve and points) and pool 
bottom coordinate versus the melting time. 0570 mm ingot, arc current J = 25 kA. 
Theory - curves, experiment - points. 



7. Conclusion. 

This article describes the mathematical simulation model of heat processes that take 
place under VAR. The discretization has been made of the non-linear heat conductivity 
equation through the Marchuck's method of integral identities and of non-linear 
boundary conditions through the final volume method. The solution algorithm in 
question has unconditional numerical stability and is applicable for tasks with non- 
commuting differentiation operators dependent on time. 

For the purposes of the simulation model testing we calculated the liquid metal 
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Figure 6. Location of the pool open surface (upper curve and points) and pool 
bottom coordinate versus the melting time. 0435 mm ingot, arc current J = 15 kA. 
Theory - curves, experiment - points. 
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pool depth under VAR melting of titanium alloy VT3-1. It was observed that the 
model simulates the liquid metal pool (depth and profiles) rather adequately during the 
whole melting process for different ingot diameters and different current strength. The 
relative error in determining the liquid metal pool profile does not exceed 15%. This 
would enable to use the model in future to calculate various melting modes in order to 
assess such parameters as: liquid metal pool depth, liquid metal pool volume, width 
of two-phase zone, temperature gradient, isotherms travel rate and local solidification 
time. 
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